Accurate effective Hamiltonians via unitary flow in Floquet space 
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We present a systematic construction of efTective Hamiltonians of periodically driven quantum 
systems. Due to an equivalence between the time dependence of a Hamiltonian and an interaction 
in its Floquet operator, flow equations, that permit to decouple interacting quantum systems, allow 
us to identify time-independent Hamiltonians for driven systems. With this approach, we explain 
the experimentally observed deviation of expected suppression of tunneling in ultra-cold atoms. 



The idea to use well-controllable quantum systems 
for simulations to explore physical phenomena has cre- 
ated big expectations to answer questions that exceed 
our computational and analytical means. In particular, 
strongly-correlated many-body states, as they typically 
occur in flat-band systems [lj-|3|, readily push us to our 
limitations of classical simulations. 

Suitably driven systems can be used to mimic the dy- 
namics of essentially any Hamiltonian, but the precise 
identification of such an effective Hamiltonian for given 
driving parameters is a big theoretical challenge. This 
can be exemplified by the dramatic increase of complex- 
ity in going from the static to the driven two-level system. 
The former is a standard textbook toy model, the latter is 
exactly solvable only in a few exceptional cases Q- The 
effective Hamiltonian of a given driven system is typi- 
cally found in an approximate manner. The deviations 
between the actual and approximated dynamics accumu- 
late in time and become significant for sufficiently long 
times. In order to perform precise quantum simulations 
it is therefore crucial to develop tools that allow one to 
systematically construct effective Hamiltonians with high 
accuracy. 

As a prominent example of driving-induced effective 
dynamics, we highlight shaken optical lattices [5|, 0, 
which permit the engineering of the tunneling of ultra- 
cold bosonic atoms confined in an optical lattice by ap- 
propriately adjusting the driving parameters. This yields 
e.g. dynamical localization [a, 0| and provides a promis- 
ingroute towards the simulation of artificial gauge fields 
[3,0]. However, despite the proven success of the usually 
considered effective Hamiltonian for high driving frequen- 
cies, substantial deviations from theoretically predicted 
dynamical localization have been observed in many-body 
samples for moderate driving frequencies [6|, which is in 
striking contrast to the single-particle case, where the 
exact dynamical localization occurs irrespectively of the 
driving frequency [13, UH ■ 

In this Letter, we introduce a new approach to derive 
effective Hamiltonians merging the concepts of Floquet 
theory [l3| and flow equations (FE) [l3|. With this, we 
provide an explanation for the experimentally observed 
deviations (a] from the theoretically predicted suppres- 
sion of tunneling [5| . 

The starting point to arrive at an effective Hamiltonian 



is a periodically driven Hamiltonian H{t) ~ H{t + T). 
After full cycles of the driving, i.e. t = nT with in- 
teger n, the time-evolution operator, U{t) =- Texp { — 

i L H{t')dt') (T denotes the time-ordering operator), 
can be written as U{nT) = ^-^H^anT ^ which defines the 
effective Hamiltonian fleff [iJl- The actual dynamics 



will follow the dynamics C/eff(t) 



-jffoi 



induced by 



the effective Hamiltonian only stroboscopically, but in 
the regime of fast driving, where w exceeds the relevant 
scales of H{t), the effective dynamics is a good approx- 
imation also for t ^ nT. The deviation between exact 
and effective dynamics defines the unitary 



UF{t)^Uetf{t)uHt) 



(1) 



Since U{t) coincides with UeS ait ~ nT, Upit) is periodic 
with period T, and equals the identity at multiples of the 
period. 

Given Eq. (jT]), the effective Hamiltonian can be 
obtained from the Schrodinger equation, i(dtU{t)) = 
H{t)U{t) {h = 1), and reads 

i^cff - UF{t)H{t)UF{t)^ - lUF{t){dtUF{t)^). (2) 

Thus, the effective Hamiltonian is found after perform- 
ing a periodic time-dependent unitary transformation, 
UF{t), such that the resulting transformed Hamiltonian 
is time- independent and C/f(0) = 1. In practice, finding 
this exact unitary transformation is an extremely diffi- 
cult task. Here we present a method, using an unconven- 
tional approach, that allows one to systematically obtain 
the effective Hamiltonian up to a required accuracy. 

We use the framework of Floquet theory [IJ, [l5| 
which asserts that the Schrodinger equation with a time- 
periodic Hamiltonian H{t) = H{t + T) has a complete 
set of solutions |0fc(i)) that decompose into a phase 
factor and a time-periodic state-vector, i.e. {(pkit)) = 
e""'^\uk{t)) with \uk{t)) = \uk{t + T)). Because of their 
periodicity, the state vectors |ufe(t)) can be expanded in a 
discrete set of periodic functions fnit), which are vectors 
in the space of time-dependent functions defined in the 
interval [0,r). In the following we will use the functions 
^muit g^g i^g^gjg g^jj(j associate with each such function a 
state vector \n) in a Hilbert space Ht. 

Any T-pcriodic operator A{t) = J^r, 



be mapped to an operator in 'Floquet space' 



A = J2^^ 



(3) 



where the Fourier components A„ = ^ /q A(i)e~*""* act 
on the Hilbcrt space of the actual system, and the (T„ act- 
ing on Ht are defined by cr„j|?i) = |?i + ?7i)[fQ|. Similarly, 
the time-derivative —idt is associated with 



V ^l(g) A ,with A = y^a;n|n)(n| 



(4) 



in the present operator notation, so that the Floquet op- 
erator K{t) ^ H{t) ~ idt is mapped to 



/C = ^ i/„ (g) o-„ -h 1 (g) w ^ n|7i) (r 



(5) 



= Hq®1+1® a + X! Hn®(Tn- (6) 

Formally, this is equivalent to a time-independent Hamil- 
tonian of a composite system with a Hamiltonian /Cq for 
the individual components and an interaction A^int- Be- 
cause of this analogy, techniques to treat interactions in 
time-independent Hamiltonians arc applicable. 

If H[t) was time-independent all Fourier components 
Hn but the static would vanish, so that the interac- 
tion term /Cint would vanish as well. Consequently, a 
static Hamiltonian is equivalent to a non-interacting sys- 
tem in the present framework. Our goal is, therefore, 
to find an operator Uc that corresponds to a periodic 
unitary transformation Uc{t) according to Eq. ([3]), such 
that the transformed Floquet operator describes two non- 
interacting systems, UcKLll = He® 1 + 1 ® A. Once 
such a transformation is found, the sought transforma- 
tion reads Upit) = U^{0)Uc{t), and the effective Hamil- 
tonian is given by iJeff = U^{0)HcUc{0). 

We target the required block-diagonalization of the 
Floquet operator with the method of flow equations 
|13l [iTf , which is considered a generalization of conven- 
tional scaling approaches and is based on a unitary flow 
that makes the Hamiltonian increasingly diagonal |18l - 
|20[ . The method defines a family of unitarily equivalent 
time-independent Hamiltonians related to each other by 
a continuous parameter I 



dHjl) 
dl 



^[V{1),H{1)], 



(7) 



where ri{l) is the anti-hermitian generator, ri{iy = —ij{l), 
of a unitary transformation. The boundary conditions 
are such that H{1 = 0) coincides with the given Hamil- 
tonian, and 'i]{l) needs to be chosen such that H{1 -^ cx)) 
is in the desired form, i.e. typically diagonal or block- 
diagonal. The canonical approach [l3| to eliminate an 
interaction i?int of a Hamiltonian H = Hq + iJjnt is to 
define the flowing Hamiltonian H{1) = Ho{l) + iJi,it(0 
and the corresponding generator ri{l) ~ [Ho{l), Hint(l)]. 



The FE method is typically used to decouple an inter- 
acting many-body system, e.g. the spin from the bosonic 
bath in the spin-boson model [21|. Here, on the other 
hand, we will use it to remove the interaction /Cint hi 
Eq. ([6]). For this purpose, we will define a flowing Flo- 
quet operator /C(/) = ICo{l) + /Cint(0 ^.nd apply Eq. ([7]) 
analogously. Additional care is however necessary in the 
choice of generator to ensure that the unitary transforma- 
tion Uc corresponds indeed to a periodic time-dependent 
transformation. This is the case exactly if U^. is invari- 
ant under the symmetry transformation S ~ 1 (g cti, 
i.e. SUcS'' = Uc, and exactly generators of the form 
Yin VniO ® cTn preserve this property. 

The generator of interest for our purposes reads 
[23, /Cint (0]- This generator will induce a flow (dynam- 
ics with the flowing parameter /) that comes to an end 
if the interaction commutes with 23. This, in turn, im- 
plies that the interaction is trivial in "Ht, J-e. it is of 
the form Hq eg) 1 and the decoupling has been achieved. 
Eq. d?]) defines a infinite set of differential equations, so 
that an exact solution can be found only in very excep- 
tional cases. We will therefore strive for a high-frequency 
expansion, where this set of equations is truncated at a 
given power in 1/w. This requires a modification of the 
generator [20| as discussed in appendix El In appendix [Bl 
we also discuss the driven two-level system for explana- 
tory purposes, and reproduce [2^ the energy shift up to 
fourth order in \/uj. Here, however, we will focus on 
the shaken optical lattice in order to address the above 
mentioned question of suppressed tunneling. 

The Hamiltonian of the one-dimensional shaken op- 
tical lattice can be written, in the co-moving reference 
frame [5|, as H{t) = Hg + Hd{t), with the Bose-Hubbard 

modelffs = Y.iJi4^i+i+4+i^i) + UJ2i'riMi-^) (with 
periodic boundary conditions) and an additional driving 
term Hd{t) = if cos(ajt) ^^ m^ that describes the shak- 
ing. J denotes the hopping matrix element (HME) be- 
tween nearest-neighbor sites and U is the on-site inter- 
action energy. The operators c^ are the usual bosonic 
annihilation (creation) operators satisfying [c2;,cj] = Sij 
and Hi = c-Ci. In lowest order in l/w we find the desired 

unitary C/i? (t) to read C/^ (t) = exp (i-^ sin(a;f) ^ jrij), 
so that the transformed Hamiltonian reads 



(0) 



Hit)^H':^>+6Hit), (8) 

with the previously known [5| effective Hamiltonian 
^cff = E.^'''(clc.+i + cl+ic.) + Uj:^n,{n, - I) and 
m{t) = E»('5+(i)c!Q+i +5"(0c1+iQ). The effective 
HME reads J°*^ = JJq{K/uj) in terms of the zeroth order 
Bessel function and 5^{t) = J(e±'-^ "'"("*> - Jq{K/uj)) 
is a small deviation. This effective Hamiltonian, H^^ , 
is a good approximation for sufficiently large driving, 
w » {J,U} 0. Experimentally, however, deviations 
from this limiting case have been observed and, as we 
shall see, a systematic improvement of _ffoff for a finite 
driving frequency w permits to explain these deviations 



very well. 
With Eq. 



(0) 



as starting point, our approach yields the 
ni+i)ci+i + h.c. , (9) 



effective Hamiltonian HeS = H^ 



including effects ~' l/i^, as described in more detail in 
appendix [C] As a qualitative change as compared to the 
lowest order effective Hamiltonian ff^g , there is tun- 
neling interaction dependent on site occupation, whose 
rate /3(X/c.) = 2^^^, J^„,_^{K / io) / {2m - 1) is given 
in terms of the m-th order Bessel functions, Jm- 

For a single particle, H^^ and all higher order terms 
vanish and our result predicts iJoff = X]j '^°*('^I'^«+i + 
CiA.iCi) as the exact effective Hamiltonian independe ntly 
of the value of J/w, in agreement with Ref. [10, lllf . 
Thus, for a single particle and even for a slow driving, the 
exact suppression of the tunneling is expected at multi- 
ples of the driving period whenever K/uj coincides with 
the zeros of Jq. With several particles, however, the ex- 
act suppression of the tunneling is not possible in general, 
since H^J and all higher order terms do not necessarily 
vanish, and only an approximate suppression in the large- 
frequency regime can be obtained. For moderate driving 
frequencies uj < 2 J [6|, H^^ is of comparable magni- 
tude as -ff^ff , so that Ec^. ^ suggests a substantially 
enhanced tunneling compared to the one predicted by 
i/^ff . Indeed, as shown in Fig.[Tl the tunneling rate pro- 
portional to (3 is close-to-maximal for driving strengths K 
at which the rate JjTo vanishes, in agreement with the ex- 
perimental observations Q. That is, in those cases W^^ 
is no longer a small higher order correction, but actu- 
ally describes the dominant tunneling mechanism. This 
shows that particular care is required in the identification 
of effective Hamiltonians and that supposedly high-order 
contributions can take over a dominant role. 

An enhancement of the tunneling rate would also ap- 
pear in a more accurate description of the optical lat- 
tice system than the Bose-Hubbard model, where next- 
nearest-neighbor tunneling is neglected assuming that 
the trapping potential is sufficiently deep [3, [2j] . H^f^ , 
however, appears also for deep lattices and, as argued 
above, its signatures can be experimentally observed even 
with a deep trapping potential, where the Bose-Hubbard 
model is an excellent approximation. 

The importance of accurate effective Hamiltonians and 
their derivation with the present method are by no means 
limited to the suppression of tunneling discussed above. 
For example, in flat band systems a small interaction 
term results in the emergence of a strongly correlated 
states [2J] , and only the correct identification of such a 
seemingly small correction to the effective Hamiltonian 
will permit the correct prediction of such states. 

Since the FE have proven very valuable in the treat- 
ment of non-perturbative effects, they bear great poten- 
tial for situations in which a clear separation of scales, as 
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FIG. 1. Rates for nearest-neighbor and tunneling interaction 
as function of the driving parameter K/u. Nearest-neighbor 
tunneling is proportional to the Bessel function Jq{K/u}) 
(solid line). The higher order process of tunneling interac- 
tion (see Eq. @) has a rate proportional to l3{K/uj) (dashed 
line). For values K/uj where Jo{K/uj) = 0, so that the tun- 
neling is expected to be suppressed, the rate for the tunneling 
interaction is close to maximal, and vice versa. Particularly 
clean realizations of Bose-Hubbard models are obtained for 
/3 = 0, as indicated with an arrow for K/lj ~ 4.3. 



assumed hitherto, is no longer valid. That is, the method 
presented here could deal with driven strongly interact- 
ing systems and identify driving schemes that simulate 
6.17. three-particle interactions. 

As a final remark, we would like to mention that 
our approach shows an equivalence between the time- 
independent and the time-dependent FE |25| . In our 
treatment, we use the framework of time- independent FE 
to treat time-dependent systems. Translating our anal- 
ysis from Floquet space back to the framework of time- 
dependent operators, one neatly reproduces the formal- 
ism of time-dependent FE. This equivalence is completely 
general and not restricted to periodic Hamiltonians, since 
one can always treat the time window of interest as the 
fundamental period of driving. That is, our approach also 
permits to translate all the existing expertise on genera- 
tors of time- independent FE to the much less mature field 
of time-dependent FE. In particular, the time-dependent 
canonical generator [25| appears as a natural extension 
of the tinre-independent canonical one. 

Financial support by the European Research Council 
is gratefully acknowledged. 



Appendix A: General approach 

Here, we will sketch how flow equations can be used to 
find transformations Upif) that yield effective Hamilto- 
nians. The starting point is the Floquet operator defined 
in Eq. (6) in the main paper 

/C = Ho ® 1 + 1 «) A + ^ i:/„ ® cr„ 

= -Ho+2? + /Cint, (Al) 



where Hq ~ Hq (g) 1. The largest energy scale of the sys- 
tem is CO, i.e. the prefactor of P, and we seek a transfor- 
mation such that the interactions in the transformed Flo- 
quet operator are of order l/w" with a predefined power 
n. In order to do so, we need a properly parametrized 
flowing Floquet operator /C(/) and generator ?7(/), which 
can be obtained in an iterative manner as described in 
the following. 

We start with an expansion of the flowing operators 



m 



(A2) 
(A3) 



i=l 



in a suitably chosen set of single system operators {Ol } 

and interaction operators {O, " }. In order to keep the 
problem tractable, this set is chosen rather small. In par- 
ticular, it is incomplete. The first ansatz for the generator 
reads 



ry(0 = [I?,/Ci„t(0] 



(A4) 



The flow equation defines differential equations for the 
parameters ai[l) and bi{l) and can be solved in lowest 
order in I/w, that is under the omission of subordinate 
terms, provided that the parametrization is chosen well. 
Essentially always, however, is a truncation necessary 
due to the expansion into an incomplete set of opera- 
tors {0\ } and {Cj }, but an inspection of the terms 
that are dropped permits to improve the parametriza- 
tion systematically: any contribution to the flow equa- 
tion that can not be spanned by {O] } and {O^ } 
defines a new operator that can be added to these sets. 
The parametrization of Ho(0 in Eq. (jA2p then takes a 
more general form (the sum includes more terms). 
A more general ansatz for the generator 



77'(0 = [P,/Ci„t(/) 



^c,({a,;(0,6,(/)})of*>i 






i-\-7n\ 



(A5) 



with suitably chosen functions q permits to solve the flow 
equation with a truncation in a higher order in the l/uj 
expansion than possible before. This process of extend- 
ing the sets {0\ '} and {Oj * } can finally be repeated 
until the desired accuracy is reached. In fact, interaction- 
like terms of n-th order can be typically neglected if an 
effective Hamiltonian in n-th order is sought, because a 
further refinement of the generator that allows to solve 
the flow equations in n-th order would result in an n-ff-st 
order correction to 7^o(0 only |20 |. 

In order to describe this procedure in more detail, we 
apply this formalism to the driven two-level system in ap- 
pendix [B] and show that it neatly reproduces the known 
high-frequency expansion j22| . In appendix[Cl we demon- 
strate that this method indeed permits to identify effec- 
tive Hamiltonians with high accuracy also for interacting 
many-body systems. 



Appendix B: Driven two-level system 

The driven two-level system is described by the Hamil- 
tonian 



H{t) = cjo-Pj + 2f^cos(a;<)P^, 



(Bl) 



where Pi, with i = x,y, z, satisfy the commutation rela- 
tion [Px, Py] ~ iPz and cyclic permutations. The Floquet 
operator, /C, associated with H{t) reads 



/C = uoPz (g) 1 + 1 ® A 
+ flPx®{ai +o--i). 



(B2) 



The operators cr„ on the second component, defined in 
Eq. (3) in the main paper, satisfy the relations (T,i(T,„ = 
CT„+m, [(T„,(T,„] = and [A, (T„] = lunon- We assume the 
regime lo ^ {^o, f^} and choose the ansatz 

/C(0 = a{l)Pz ® 1 + 1 ® A + b{l)Px ® (cti + cr_iiiB3) 

with the boundary conditions a(0) = uq and 6(0) = O. 
There are only two flowing variables a{l) and 6(/), that 

is, the sets {0\ } and {0^ } are comprised only of the 
single elements Pz®l and Px ® (ui + a-i) respectively. 
With the generator 

771 (0 = [P,/Cint(0] = '^b{l)Px ® (ai - a_i), (B4) 
the flow equation reads 
I d/C(0 



w^ 



dl 



-b{l)Px®{(yi+(r-i) 



-a{l)b{l)Py ® (cti 



^-i 



(B5) 



Eq. (|B5|1 can only be solved, if the first order term 
— a{l)b{l)Py (E) ((Ti — (T_i) is neglected, since it can not be 
expanded in term of the operators considered hitherto. 
In this lowest order, one obtains the differential equation 
-^ — —uP'bil), which yields an exponential decay of b 
with I, such that the interaction vanishes for / — >■ oo as 
expected. Nevertheless, as argued above, one will obtain 
a better approximation to the exact flow equations, if the 

term Py ® (cti — cr„i) is added to the set {Of" }. 
With the modified generator 



772(0 = ujb{l)Px (E) {ai -cr_i) 
+ ci{l)Py®{ai +cr_i) , 

the flow equation becomes 

I dlC{l) 



(B6) 



dl 



= -b(l)Px®{ai+<7^i) 



1 



{ia{l)b{l) + ci{l))Py (g) (cti - cr_i) 



-2^ci{l)b{l)P,®l 

+ -!^Ci{l)a{l)Px®{cTi+a^i) 



i 



Ci{l)b{l)Pz® {(72 +(7-2). 



(B7) 



With the choice Ci{l) = —ia{l)b{l) all the first order (^ 
— ) interaction terms cancel, but there is the second order 

a; '' ' 

interaction term Pz(8)(cr2+fT_2)- Since it is off-diagonal, it 
yields no contribution to the diagonal part of the Floquet 
operator in order l/w^ and it will be eliminated by a 
suitable choice of 77 in the order l/w^. 

A second order effective Hamiltonian is therefore given 
in terms of the truncated flow equation 



1 rf/C(0 



aHl)^ 






w^ 



(B8) 



The solutions with boundary condition /C(0) = /C satisfy 



^2 

2 . 



a{l ^ 00) = LUo{l :^)+0{n'^/uj^) and (B9) 

b{l^ 00) ^0 + O{n^/Lj^) (BIO) 

and lead to the constant Hamiltonian 






(Bll) 



Solving the flow equations in fourth order requires 
adding the operators P^ ® (0-2 — o'-2), Py ® (0-3 + cr-3) 
and Py ® (ci + CT_i) to the set {O^ " } and taking the 
generator 

774(0 = ujb{l)P^ (g) (cTi - cr_i) - ia{l)b{l)Py (g) (crj + cr_i) 
a{l)bm) 



2uj 
.a{l)b^{l) 



Pz ® (ct2 - (y-2) 



Py ® (era + Cr_3) 

2 Py®{(Tl+a-l) 



6C.2 ^y 

a{l)b^{l) 

2uj 



(B12) 



where suitable choices for the functions c,; have alreay 
been made. With this ansatz, the flow equation reads 

- -b{l){\ ^ P^«)(cri+a_i) 



dl 



- [^a{l)b\l) + 



a(0&^(0 



P^® 1 



/(l/^'), 



(B13) 



where /(1/w^) are interaction terms ^ l/oj* that can be 
neglected, similarly to above. The solutions of Eq. (|B13p 
with the proper boundary conditions satisfy 






&(Z^oo) =0 + O(OVw'^) 



— ) + o(f^V-^), 

(B14) 



Again, the interaction vanishes for Z — > 00, and Eq. (|B14p 
reproduces the effective level shift of a driven two-level 
system in fourth order |22l |. 

Appendix C: Shaken optical lattices 



The Floquet operator of the shaken optical lattice 
Hamiltonian reads 



^ = ^^aJ{4c^+l+4+lC■^ + UY,n^{n, - 1)) «) 1 

(CI) 



1(8) A 



where 5 yt), defined in Eq. (8) in the main pa- 
per, has been expressed as a Fourier series 5 {i) = 



JTZ=-^J±^n{K/u:)e 



JJq{K/uj) with the Tn-th 



order Bessel function. As above, we consider the regime 
of fast driving, {J, U} <C lo. We start with the following 
parametrization of the fiowing Floquet operator: 



/Ci(0 = (^J^J{4c^+l+c\^^c,) + U^n,{n, - 1)) 8) 1 

i i 

-t- 1(8) A 

+ X! ^™(0'^(^m4<=»+l +^-™4+lCi) (XlCTm. (C2) 
2,7717^0 



According to our discussion above, the generator should 
read 



77i = [2?,/Cint(0] (C3) 

UJmb,n{l)J{Jmc\ci+l + J-mc\^iCi) ® dm- 



E 



(C4) 



With this choice, however, one would obtain Upit = 0) 7^ 
1. In order to avoid the additional step of introducing 
a constant Hamiltonian rather than identifying i/cft di- 
rectly, we choose a slightly modified generator 

'71=^ WTO^m(0'^(^m4c7+l+^-mC- + lCj)®(o-m~l)- 

(C5) 
With this choice, the flow equation becomes 



1 u/Clf/) rt^ ^ X~^ n\ I rr \f \ rr \ I \ \ 

— — -j^ — = 2 2^ Tnbm(l)[JmCl (ni - n^+ijCi+i + J-mC,-_^i(nj+i - ni)Ci) ® 1 



a;2 d; 



■J.mT^O 



UJ 



- 2 ^ mbm{l){Jmc\{ni - nj+i)ci+i + J-mcJ_|_i(ni+i - ni)ci) ® cr„ 



i.TnT^O 






(C6) 



In contrast to the case of the driven two-level system above, it is necessary to also expand the set [O] }. Specifically, 
the tunneling interaction term needs to be added, and will appear in the resulting effective Hamiltonian. With the 
extended set {Ol } the flowing operator reads 

/C2(/) = 1 A + {^JoJ{c\c,+^ + c\^^c,) + U^n,{n, - 1)) «) 1 

i i 

+ 2^ {a+{l)c\{ni - nj+i)ci+i + a_(/)c.|^^(ni+i - ni)c^) ® 1 






(C7) 



Again, the first order interaction terms in Eq. (jC6p can be neglected, so that there is no need to introduce a modified 
generator, but one directly obtains the truncated flow equation 



— — J. — = 2 ^ mbm{l){jmcl(ni - nj+i)cj+i + J'-mcJ_,_i(ni+i - n,;)ci) ® 1 



uj^ dl 



i;m^{) 



- J ^ m'^bm{l){Jr,ic\ci+l + J-mc\^iCi) ® Um- 

The solutions of Eq. (jCSp with the initial conditions a±(0) = satisfy 

a±{l ^ cx)) = ±—l3 + OiJ^uj^) and 
a; 



bmil^Oo) =0 + O(J7w^) , 



(C8) 



(C9) 
(CIO) 



with/3 = E„.o^^-^ = 2E,^=i 



OO J2rr^-l{K/u:) 



Consequently, IC2 satisfies 



UJ- 



IC2{1 -!■ cx)) = ('^JoJicjc^+i + cj^iCi) + C/^?ij(ni - I) + 2(3 ^c|("-» ^ nj+i)ci+i + h.c.j ® 1 



+ 1®A + 0(JVw^) , 

which results in the effective Hamiltonian given in Ec^. ^ in the main paper. 
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